c=[-9;-8;-50;-19];
A=[3 2 3 4;0 0 2 2.5];
b=[18;3];
lb=[0;0;0;0];
ub=[Inf;Inf;Inf;Inf];
[x,fval]=linprog(c,A,b,[ ],[ ],lb,ub)


Max = 9*x1 + 8*x2 + 50*x3 + 19*x4;
3*x1 + 2*x2 + 10*x3 + 4*x4 <= 18;
2*x3 + 2.5*x4 <= 3;


Max = 9*x1 + 8*x2 + 50*x3 + 19*x4 + x5;
3*x1 + 2*x2 + 10*x3 + 4*x4 + 3*x5 <= 18;
2*x3 + 2.5*x4 + x5 <= 3;
